SNSN-323-63 
July 19, 2012 



Non Perturbative One Gluon Exchange Potential from 
Dyson-Schwinger Equations 



P. Gonzalez 1 , V. Vento 1 

1 Departamento de Fisica Teorica -IFIC. Universidad de Valencia- C SIC. E-^6100 

Burjassot (Valencia), Spain. 

V. Mathieu 2 

2 ECT* , European Center for Theoretical Studies in Nuclear Physics and Related 
Areas, Strada delle Tabarelle 286, 1-38123 Villazzano (TN), Italy 



Recent progress in the solution of Dyson-Schwinger equations of QCD 
allows for a non perturbative evaluation of the One Gluon Exchange 
(OGE) interaction. We calculate the interquark static potential for heavy 
mesons by assuming that it is given by this OGE interaction and we apply 
it to the description of charmonium. 
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1 Introduction 



The development of non-perturbative techniques in Quantum Chromodynamics (QCD) 
is starting to allow for a description of the hadron spectrum from first principles. 
Among these techniques Lattice gauge theory (TJ [2] constitutes a non-perturbative 
regularisation scheme which may provide accurate numerical solutions of the the- 
ory. The accuracy of lattice results has been tremendously improved during the 
past decade with the availability of more powerful computers [31 HJ. Therefore, lat- 
tice results are considered in many instances the data which other non perturbative 
approximations try to reproduce. One of these alternatives is the approximate resolu- 
tion of the Dyson-Schwinger Equations (DSE) of QCD, a non-perturbative approach 
which has progressed considerably in the last ten years. This approach, more analyt- 
ical, has led to a very appealing physical picture establishing that the QCD running 
coupling (effective color charge) freezes in the deep infrared. This property can be 
best understood from the point of view of a dynamical gluon mass generation that is 
a purely non-perturbative effect [5j [6j [T] . 

The aim of this presentation is to investigate, following reference [8], the form of 
the OGE static potential from DSE and compare it to the static potentials derived 
from lattice calculations. The application of these potentials to the description of 
quarkonia will be discussed. 

2 One Gluon Exchange Potential from Dyson 
Schwinger Equations 

It is well established by now that the QCD running coupling (effective charge) freezes 
in the deep infrared what may be understood in terms of a dynamical gluon mass. 
Even though the gluon is massless at the level of the fundamental QCD Lagrangian, 
and remains massless to all order in perturbation theory, the non-perturbative QCD 
dynamics generate an effective, momentum- dependent mass, without affecting the 
local SU(3) C invariance, which remains intact [5j [TJ [9] . At the level of the Dyson- 
Schwinger equations (DSE), solved by using a PT (pinch technique) - BFM (back- 
ground field method) truncation scheme in the quenched approximation (no quark 
loops), the generation of such a mass is associated with the existence of infrared fi- 
nite solutions for the gluon propagator. Such solutions may be fitted by "massive" 
euclidean propagators of the form A _1 (g 2 ) = q 2 + m 2 (q 2 ) where m 2 (q 2 ) depends 
non-trivially on the momentum transfer q 2 . 

One physically motivated possibility is the so called logarithmic mass running, 
which is defined by 
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where A is the QCD scale and mo, the gluon mass at q 2 — > 0, as well as p and 6 are 
constants to be fitted. For this fitting the calculated DSE propagator is compared to 
the one obtained from lattice. 

On the other hand the non-perturbative generalization of ot{q 2 ), the QCD running 
coupling, comes in the form 
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where /3o = H — 2ri//3 being n/ the number of active quark flavours. 

Note that its zero gluon mass limit leads to the LO perturbative coupling con- 
stant momentum dependence. The m(q 2 ) in the argument of the logarithm tames 
the Landau pole, and a(q 2 ) freezes at a finite value in the IR, namely a(0) = 

4tt [/3 ln(^)] _1 . 

From these expressions for the propagator and the coupling we can derive the One 
Gluon Exchange interaction between static charges, see FigHJ 




Figure 1: One gluon exchange interaction. 

The OGE static potential is related to the Fourier transform of the time-time 
component of the full gluon propagator in the following way 

V (r) --c.(— 4?m(k2) c *» - 32nCF f°dk k a(k2) sinfkr) (3) 

where Cf is the Casimir eigenvalue of the fundamental representation of SU (3) \Cp = 
(N 2 — 1)/2N for SU(N)\ and the bold terms, k and r stand for 3-vectors. 
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In Fig. |2] we show the finite part (up to a constant) of the potential, Eq[3j derived 
from the DSE with the definitions in Eqs [2] and HJ We have chosen the following 
range of parameters: m ~ 360 - 480 MeV, p = 1. — 4., <5 = l./ll EH Q3] which 
provide a good fit to the lattice propagator. The value of A here has been taken to be 
300 MeV. In order to adjust the behavior at the origin to the data we have used /?o 
corresponding to rif = 4 flavors. To do this appropriately one should introduce the 
running of the quark masses, which are at present not well known. However, since 
asymptotically the masses run to zero, our way of proceeding achieves the correct 
(perturbative) strength of the potential at low r. The potential describes well the low 
radial behavior, by construction, flattens at large r going asymptotically to zero and 
never becomes positive. 

For comparison we have used i) a Cornell type potential (whose form is derived 
from quenched lattice calculations [121 03] ) 



containing the perturbative expectation plus an additional linear term and ii) a 
screened type potential (whose form was derived long time ago from an unquenched 
lattice calculation [14] ) 



The values used for the parameters are based on spectroscopy (see next Section). 

We should realize that the additive infinite self-energy contribution associated 
with the static sources should be removed from the calculated OGE potential [12] . 
In lattice QCD this is done normalizing the potential such that V(tq) = where 
is the Sommer scale [15] with a typical value around 0.5 fm. We proceed the same 
way but for phenomenological purposes we take the substraction point at the zero 
point of the potentials used in Table 1, which happens to be at r ~ 0.35/m for the 
parametrizations we are using. The result of this procedure is shown in Fig. |3j This 
procedure increases the value of the potential without changing its shape. 

As can be checked our OGE potential resembles the screened and not the Cornell 
potential. We have shown that none of the parameters plays a fundamental role in 
the structure of the DSE OGE potential, since the structure does not vary when we 
vary them. All values discussed lead to the same qualitative features for the potential. 
It is evident that there is no way to reproduce the large r behavior of the Cornell 
potential by changing the parameters. The DSE OGE potential flattens and becomes 
basically constant similar to the screened potential. If we assume that the OGE 
interaction is the main source of the dynamics, this result is quite surprising since the 
approximations used to find the solution to the DSE do not contain quark loops and 
therefore they incorporate no mechanism for screening, i.e. some mechanism derived 
from the breaking of the string [T4| [16] . 
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Figure 2: The DSE OGE potential, i.e. the Fourier transform of the massive One 
Gluon Exchange, is plotted for the range m ~ 360 — 480 MeV and p ~ 1 — 4 with 
rif = 4 (Pq = 25/3). For comparison we plot the Cornell and screened potentials with 
parameters a = 0.52, y/b = 427 MeV and 7 = 0.38 fin -1 . 

Therefore we should conclude that either some correction to the calculated OGE 
is lacking (this could come from the use of a higher order truncation scheme for the 
resolution of the DSE, from vertex corrections to the coupling ...) or that the incor- 
poration of multigluon effects is essential for deriving a consistent large r behavior. 
In this respect an ad hoc explanation based on OGE dominance plus additional cou- 
pling corrections has been recently proposed along the following arguments jTT] . Let 
us consider that there were a nonperturbative vertex correction to the strong effective 
charge in the form 

a CO nf{q ) = 
with c a constant to be fitted, so that 

atotai{q 2 ) = a conf (q 2 ) + a(q 2 ) 

where the last term on the right hand side would be given by Eq[2j Then it would 
be possible to fit not only the lattice propagator but also the Cornell type lattice 
potential and a fully consistent scheme in the quenched approximation would arise. 
To go further one should incorporate quark loop effects. According to reference [TS] 
the former vertex correction could be effectively modified. If this modification were 
parametrized through a cutof (s) in the form 

rfA 4 
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Figure 3: The DSE OGE potential after the Sommer substraction, with the same 
parameters as in the previous figure, is plotted. For comparison we plot the Cornell 
and screened potentials. 



being d a constant, the resulting DSE OGE potential would be like the screened 
potential. In this manner an explanation of our previous results compatible with 
one gluon exchange dynamical dominance could come out. Next we shall adopt this 
interpretation in the phenomenological application of the SDSE OGE potential to 
charmonium. 



3 Quarkonia description 

For sufficiently heavy quarks, one may hope that the bound state problem becomes es- 
sentially non-relativistic the dynamics being controlled approximately by a Schrodinger 
equation with a static potential. It should be emphasized that the derived potentials 
do not contain spin-dependent terms what make them more reliable when these terms 
do not play a major role. This is expected to occur for high excited (large-sized) 
states since spin corrections are short-ranged. For the low lying states data show 
that spin-spin splittings between spin triplet and spin singlet states should be rele- 
vant {mj/^ - m Vc = 117 MeV, m^ 2s - m Vc ^s) = 49 MeV, m T (i s ) - m Vb = 69 MeV, 
...). By considering that the perturbative spin-spin correction is in absolute value 
three times bigger for singlets than for triplets, the radial potential approach could 
be taken as an approximate description of spin triplet states. Other spin-dependent 
corrections (spin-orbit, tensor) may be playing some role. It is worth to point out 
that perturbative spin-orbit and tensor splittings cancel in the centroids of p waves 
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which consequently may be used as "data" for comparison with the radial approach 
results. Further relativistic corrections are expected to be more important for char- 
monium than for bottomonium. Therefore the radial approaches are better suited 
for the study of bottomonium. One should not forget though that in the applica- 
tion to quarkonia the parameters entering the expressions of the potentials have an 
effective character since their values may be implicitly incorporating non considered 
corrections. 

The charmonium spectrum with well established quantum numbers, correspond- 
ing mostly to J PC = 1 resonances produced through ISR (Initial State Radiation) 
processes, is analized in Table 1. The calculated spectrum from the DSE OGE poten- 
tial is compared to experimental data and to a Cornell like potential calculation (the 
choice of charmonium instead of bottomonium makes clearer the differences between 
both calculations). The parameters for the Cornell potential have been chosen within 
the conventional spectroscopic range a ~ 0.51-0.52 and Vb ~ 412-427 MeV [201 EE] 
(note that concerning the description of quarkonia masses the addditive constant in 
the potential can be absorbed in a renormalization of the quark mass). 



n r L 


Mcornell 


M DSE 


M PDG 




MeV 


MeV 


MeV 


Is 


3069 


3151 


3096.916 ±0.011 


2s 


3688 


3660 


3686.09 ± 0.04 


Id 


3806 


3761 


3772.92 ±0.35 


3s 


4147 


4004 


4039 ± 1 


2d 


4228 


4070 


4153 ±3 


As 


4539 


4273 


4263t^ 


3d 


4601 


4321 


4361 ± 9 ± 9 


5s 


4829 


4487 


4421 ± 4 


Ad 


4879 


4526 




6s 


5218 


4651 


4664 ± 11 ± 5 


lp 


3502 


3515 


3525.3 ±0.2 


2p 


3983 


3886 





Table 1: Calculated masses, Mcomeii and Mdse, from the Cornell and DSE OGE 
potentials. For the Cornell potential a = 0.52, \/b = 412 MeV and m c = 1350 MeV. 
For the DSE potential m = 345.7 MeV, p = 1 and m c = 1400 MeV. Masses for 
experimental candidates, Mp^c, have been taken from [19]. For p waves we quote 
the centroid of npo, np\ and np2 states. 

As can be checked the main difference between the two models refers to the de- 
scription of the higher excited states. The DSE model allows in the overall for a 
reasonable one to one assignment of calculated states to data (within 60 MeV differ- 
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ence) whereas the Cornell model, providing a good fit for the lower states (at most 
30 MeV difference with data), can not accomodate all the known higher energy res- 
onances but only some of them. For instance ■0(4040), 0(4160) and 0(4415) may be 
assigned to the Cornell 3s, 2d and 4s states respectively. Then other two resonances, 
cataloged in the Particle Data Group Review [19] as X(4260) and X(4360), can not 
be reproduced (this has motivated alternative, non cc, explanations for these states 
even though their properties might be understood as corresponding to cc states, see 
[22] and references therein). 

Although more complete analyses are needed before extracting any definite con- 
clusion these results seem to point out that the nonperturbative OGE potential might 
provide a well founded approach to heavy meson spectroscopy. 

4 Summary 

We have calculated the OGE static potential from an approximate solution of the 
quenched Dyson-Schwinger equations for the gluon propagator. The low r behavior 
is determined by the well know asymptotic behavior. The large r behavior is certainly 
non perturbative. The Sommer procedure, to avoid self energy effects of the static 
charges, leads to a potential which is not negative everywhere. The DSE with the 
Sommer normalization is quite similar to a screened potential form derived from 
unquenched lattice calculations. 

This suggests that non considered nonperturbative vertex corrections to the strong 
effective coupling could be responsible for the linear confining. These corrections 
might be canceled in the unquenched DSE solution providing an ad hoc explanation 
of the OGE results obtained. 

Taking for granted this explanation and assuming that the OGE interaction is the 
main source of the dynamics we have proceeded to a calculation of the charmonium 
spectrum. We have found a one to one correspondence between the calculated states 
and the experimental resonances. This makes us tentatively conclude that the nonper- 
turbative OGE interaction may provide a significant improvement in the description 
of heavy quarkonia as compared to conventional potentials based on confinement plus 
perturbative OGE terms. 
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